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Abstract 

In this expository article we collect the integer sequences that count several different 
types of matrices over finite fields and provide references to the Online Encyclopedia of 
Integer Sequences (OEIS). Section 1 contains the sequences, their generating functions, 
and examples. Section 2 contains the proofs of the formulas for the coefficients and 
the generating functions of those sequences if the proofs are not easily available in 
the literature. The cycle index for matrices is an essential ingredient in most of the 
derivations. 

2000 Mathematics Subject Classification: Primary 05A15, 15A33; Secondary 05A10, 
05A30, 11B65, 11B73. 

Keywords: Integer sequences, Matrices over finite fields, Cycle index for matrices. 

Notation 

q a prime power 

F g field with q elements 

GL n (g) group ofnxn invertible matrices over F q 

M n (q) algebra ofnxn matrices over F q 

7n(<z) order of the group GL n (g) 

7 n order of GL n (g) with q understood 

Vd number of irreducible monic polynomials of degree d over F„ 
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Sequences 


Section 


OEIS 


All matrices 


1.1 


A002416 


Invertible matrices 


1.2 


A002884 


Subspaces, g-binomial coefficients 


1.3 


A022166-188 






A006116-122 






A015195-217 


Splittings (direct sum decompositions) 


1.4 




g-Stirling numbers, g-Bell numbers 


1.4 




Flags of subspaces 


1.5 


A005329, A069777 


Linear binary codes 


1.6 


A022166, A076831 


Matrices by rank 


1.7 




Linear derangements 


1.8, 2.2 


A002820 


Projective derangements 


1.8, 2.2 




Diagonalizable matrices 


1.9, 2.3 




Projections 


1.10, 2.4 


A053846 


Solutions of A k = I 


1.11, 2.5 


A053718, 722, 725 






A053770-777 






A053846-849 






A053851-857 






A053859-863 


Nilpotent matrices 


1.12 


A053763 


Cyclic (regular) matrices 


1.13, 2.6 




Semi-simple matrices 


1.14, 2.7 




Separable matrices 


1.15, 2.8 




Conjugacy classes 


1.16, 2.9 


A070933, A082877 



1 The Sequences 

References to the OEIS are accurate as of February 1, 2006. Sequences mentioned in this 
article may have been added since then and entries in the OEIS may have been modified. 

1.1 n x n matrices over F q 

Over the field F q the number ofnxn matrices is q n . For q = 2 this sequence is A002416 
of the OEIS, indexed from n — 0. The terms for < n < 4 are 

1,2,16,512,65536. 

1.2 Invertible matrices 

The number of invertible n x n matrices is given by 

7 „(?) := \GK(q)\ = (<f - 1)(? B - ?)(?" - q 2 ) • • • (<f - q^ 1 )- 



2 



We will use 7„ in place of 7„(g) unless there is a need to be explicit about the base field. It 
is convenient to define 70 to be 1. For q = 2 the sequence is A002884, and from 7 to 75 the 
terms are 

1,1,6,168,20160,9999360. 

A second formula for 7„ is 

7» = (9 - 1) B <A4,!, 

where [n] q \ := [n] q [n — l] q • ■ • [2] q [l] q and [i] q : = 1 + q + q 2 + ■ ■ ■ + g l_1 are the g-analogs of n\ 
and i. Also, [i] q is the number of points in the projective space of dimension i — 1 over F q . 
A third formula is 

7n = (-l) n ?^(?;?)n, 

where (a; q) n is defined for n > by 

n-l 

Construct a random n x n matrix over F g by choosing the entries independently and 
uniformly from F q . Then r y n /q n is the probability that the matrix is invertible, and this 
probability has a limit as n — > 00 

For q = 2 the limit is 0.28878 . . .. For g = 3 the limit is 0.56012 .... As g — > 00 the probability 
of being invertible goes to 1. 

1.3 Subspaces 

The number of /c-dimensional subspaces of a vector space of dimension n over F q is given by 
the Gaussian binomial coefficient (also called the g-binomial coefficient) 

fn\ _ (g n - l)(g n - q)(q n - q 2 ) ■ ■ ■ (q n - q^ 1 ) 
\k) q ~ (q k -l)(q k -q)(q k -<f)---(q k -q k - i y 

Also, we define (°) = I. Other useful expressions for the Gaussian binomial coefficients are 

(A W 

\k) q [k] q \[n-k] q V 

which shows the g-analog nature of the Gaussian coefficients, and 

fn\ 7» 
\k) q lkln-kq k{n ~ k) ' 

which comes from the transitive action of GL n (g) on the set of subspaces of dimension k, 
the denominator being the order of the subgroup stabilizing one of the subspaces. 
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Entries A022166-A022188 are the triangles of Gaussian binomial coefficients for q = 2 to 
q = 24. Note that when q is not a prime power the formula does not count subspaces. For 
q = 2 the rows from n = to 6 are 

1 

1 1 

1 3 1 

17 7 1 

1 15 35 15 1 

1 31 155 155 31 1 

1 63 651 1395 651 63 1 

For the Gaussian binomial coefficients we have the g-binomial theorem [2] 

n a+**>= e « (!) 

l<i<n 0<fc<n ^ ' 1 

Summing over k from to n for a fixed n gives the total number of subspaces. Sequences 
A006116-A006122 and A015195-A015217 correspond to the values q = 2, . . . , 8 and q = 
9,..., 24. (Same warning applies to those q that are not prime powers.) For q = 2 and 
n — 0, ... ,8 the sequence begins 

1, 2, 5, 16, 67, 374, 2825, 29212, 417199. 



1.4 Splittings 

Let be the number of direct sum splittings of an n-dimensional vector space into k 
non-trivial subspaces without regard to the order among the subspaces. These numbers are 
g-analogs of the Stirling numbers of the second kind, which count the number of partitions 
of an n-set into k non-empty subsets, and the notation follows that of Knuth for the Stirling 
numbers. Then it is easy to see that 




It is not difficult to verify that the triangle satisfies the two variable generating function 
identity 

i+££{r} ^-(^v 

n>l k=l k J 1 ,n \ r>l ,r / 

For q = 2 the entries for 1 < n < 6 and 1 < k < n are 
1 

1 3 

1 28 28 

1 400 1680 840 

1 10416 168640 277760 83328 

1 525792 36053248 159989760 139991040 27998208 
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Let b n be the total number of non-trivial splittings of an n- dimensional vector space: 



n 



n 



The b n are analogs of the Bell numbers counting the number of partitions of finite sets. Then 
we have the formula for the generating function for the b n , 




For q = 2 the values of b n for n = 1 , . . . , 6 are 

1, 4, 57, 2921, 540145, 364558049. 



Bender and Goldman pQ first wrote down the generating function for the g-Bell numbers, 
which, along with the g-Stirling numbers, can be put into the context of g-exponential families 

m 



1.5 Flags 

A flag of length k in a vector space V is an increasing sequence of subspaces 

{0} = V c V 1 C V 2 C • ■ ■ V k = V. 

If dim V = n, then a complete flag is a flag of length n. Necessarily, dim V^ — i. 
The number of complete flags of an n-dimensional vector space over F q is 

l) ("l 1 ) ---(J) ^[nUn-ll.-'-Ill^Iny. 

The reason is that having selected Vi, . . . , Vi, the number of choices for Vi + i is Alter- 
natively, the general linear group acts transitively on the set of complete flags. The stabilizer 
subgroup of the standard flag is the subgroup of upper triangular matrices, whose order is 
(q — l) n g( 2 ). Hence, the number of complete flags is 

In r 1 1 

— pr = W- 

(q - l) n qW 

Define [0] 9 ! = 1. The sequence [n] q \ for q = 2 is A005329. The terms for n = 0, 1, 2, . . . , 8 

are 

1, 1, 3, 21, 315, 9765, 615195, 78129765, 19923090075. 

By looking up the phrase "g-factorial numbers" in the OEIS one can find the sequences 
[n] q \ for q < 13. The triangle of g-factorial numbers is A069777. 
The g-multinomial coefficient 



n 



K n 1 n 2 ■••n k J q [ni] g ![n 2 ] g ! • • • [n k ] q \ 

where n = n\ + ■ ■ ■ + n k , is the number of flags of length k in an n-dimensional space such 
that dim Vi = n\ + • — h rij. 
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1.6 Linear binary codes 



An [n, k] linear binary code is a fc-dimensional subspace of the space of F£ . Thus, the number 
of [n, k] linear binary codes is the Gaussian binomial coefficient (^) from section 1.3 and is 
given by the triangle A022166. 

Two linear binary codes are equivalent (or isometric) if there is a permutation matrix 
P mapping one subspace to the other. The number of equivalence classes of [n, k] codes 
is given by the triangle A076831. The early entries are identical with the corresponding 
binomial coefficients, but that does not hold in general. The rows for n = to 6 are 

1 

1 1 

1 2 1 
13 3 1 
1 4 6 4 1 
1 5 10 10 5 1 
1 6 16 22 16 6 1 

For linear codes over other finite fields, the notion of equivalence uses matrices P having 
exactly one non-zero entry in each row and column. These matrices, as linear maps, preserve 
the Hamming distance between vectors. They form a group isomorphic to the wreath product 
of the multiplicative group of F q with S n . 

There are more than a dozen sequences and triangles associated to linear binary codes 
in the OEIS. They are listed in the short index of the OEIS under "Codes." 



1.7 Matrices by rank 

The number ofmxn matrices of rank k over F q is 



"Ij (? B -l)(? B -?)-"(? B -9*- 1 ) 



k 



(q m -l)(q m -q)---(q m -q k 



i 

= (q m - l)(q m - q) ■ ■ ■ (q m - g^ 1 ) (q n - l)(g ra - g) • • • (q n - q k ~ l ) 
(q k — l)(q k — q) • • • (q h — q k ~ v ) 

(To justify the first line of the formula note that the number of /c-dimensional subspaces of F^ 1 
to serve as the column space of a rank k matrix is (^) ■ Identify the column space with the 

image of the associated linear map from F™ to F™. There are (q n — l)(q n — q) ■ ■ ■ (q n — g^" 1 ) 
surjective linear maps from F™ to that /c-dimensional image. The second line follows by 
transposing.) 

Define the triangle r(n, k) to be the number of n x n matrices of rank k over F q and 
r(0,0) := 1. Thus, 

((q n - l)(q n - q) ■ ■ ■ (q n - q k ~ r )) 2 



r(n, k) = 



(q k — l)(q k — q) ■ ■ ■ (q k — q k x ) 
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For q = 2 the entries from r(0, 0) to r(5, 5) are 



1 

1 1 

1 9 6 

1 49 294 168 

1 225 7350 37800 20160 

1 961 144150 4036200 19373760 9999360 



1.8 Linear and projective derangements 

A matrix is a linear derangement if it is invertible and does not fix any non-zero vector. 
Such a matrix is characterized as not having or 1 as an eigenvalue. Let e n be the number 
of linear derangements and define eo = 1. Then e n satisfies the recursion 

e n = e n _M n ~ l)?"" 1 + (-1) V {n " 1)/2 . 

For q = 2, the sequence is A002820 (with offset 2) and the first few terms beginning with e 
are 1, 0, 2, 48, 5824, 2887680, The sequence can be obtained from the generating function 

n>l ln r>l V q / 

The proof of this is in section 2. The asymptotic probability that an invertible matrix is a 
linear derangement is 




The asymptotic probability that a square matrix is a linear derangement is 




A matrix over F q is a projective derangement if the induced map on projective space 
has no fixed points. Equivalently, this means that the matrix has no eigenvalues in F q . If 
d n is the number of such n x n matrices, then 

n>l ln r>l V q J 

For q = 2 the notions of projective and linear derangement are the same, and the corre- 
sponding sequence is given above. For q = 3 the sequence has initial terms from n — 1 to 
5 

0, 18, 3456, 7619508, 149200289280. 
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Since two matrices that differ by a scalar multiple induce the same map on projective space, 
the number of maps that are projective derangements is d n /(q — 1). The asymptotic proba- 
bility that a random invertible matrix is a projective derangement is the limit 



,. dr. 
hm — 



"- X) In " V q r 



9-1 



r>l 

The asymptotic probability that a random n x n matrix is a projective derangement is 



lim 



n— >oo qr" \ Q 



As g — ■> oo, this probability goes to 1/e, the same value as the asymptotic probability that 
a random permutation is a derangement. 

1.9 Diagonalizable matrices 

In this section let d n be the number of diagonalizable n x n matrices over F q . Then 

n>l 7n \m>0 7 "\ 

It follows that 

d n = V • 

r • • • '"V 

„ 1+ ...+n,=n ' ni ' n <* 

For g = 2 this simplifies to 



n 



d n = J2 

The sequence for d\ to d 8 is 



In 



2, 8, 58, 802, 20834, 1051586, 102233986, 196144424834. 
For q = 3 the five initial terms are 

3, 39, 2109, 417153, 346720179. 
For arbitrary q, we can easily find d 2 , 

In 



„ 1 + ...+„ q =2 7ni " ' 7n « 

7o72 V 2 / 7i7i 
g 4 - g 2 + 2g 
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In 



1.10 Projections 

A projection is a matrix P such that P 2 = P. Let p n be the number ofnxn projections. 
Then 

l + yV_n u n = 

n>l 7n 

It follows that 

Pn= 

0<i<n Wn-i 

In the sum the term ^n/liln-i is the number of projections of rank %. Projections are also 
characterized as diagonalizable matrices having eigenvalues or 1. Thus, for q = 2 the 
diagonalizable matrices are precisely the projections, and so we get the same sequence as in 
section 1.9. From p 1 to p 8 the sequence is 

2, 8, 58, 802, 20834, 1051586, 102233986, 196144424834. 

For q = 3 the map P i— > P + / gives a bijection from the set of projections to the set of 
diagonalizable matrices with eigenvalues 1 or 2. Such matrices are precisely the solutions of 
X 2 = I. The sequence is given in the OEIS by A053846. From p 1 to p 7 the sequence is 

1, 2, 14, 236, 12692, 1783784, 811523288, 995733306992. 

There is a bijection between the set ofnxn projections and the direct sum splittings 
F™ = V © W. The projection P corresponds to the splitting ImP © KerP. Note that V © W 
is regarded as different from W © V. Hence, 



Pn = 2 + 2 



because is the number of splittings into two proper subspaces, with V © W regarded 
as the same as W © V. 
From §1.3 we have 



kj q lkln-kq k{n ky 



showing the relationship between the number of subspaces of dimension k and the number 
of projections of rank k. From it we see that t there are q k( - n ~ k ) complementary subspaces 
for a fixed subspace of dimension k. 



1.11 Solutions of A k = I 

Let a n be the number of n x n matrices A satisfying A 2 = I . Such matrices correspond to 
group homomorphisms from the cyclic group of order 2 to GL n (q). In characteristic other 
than two, the generating function for the a n is 




9 



and so 



n 




In 



i=0 



TiTn— i 



These matrices are those which are diagonalizable and have only 1 and —1 for eigenvalues. 
If we fix two distinct elements Ai and A2 in the base field F q , then a n also gives the number 
of diagonalizable matrices having eigenvalues Ai and A2. Taking the eigenvalues to be and 
1 gives the set of projections, and so we see that the number of projections is the same as the 
number of solutions to A 2 = I, when q is not a power of 2. The proof in §2.4 for the number 
of projections is essentially the proof for the more general case. The sequence A053846 gives 
the number of solutions over F 3 . 

In characteristic two, A 2 = I does not imply that A is diagonalizable and the previous 
formula does not hold. We have the formula 



to be proved in §2.5. Because A 2 = I is equivalent to (A + I) 2 = 0, the formula for a n also 
counts the number of n x n nilpotent matrices N such that iV 2 = 0. For q = 2 the sequence 
is A053722. Although there appears to be an error in the formula given in the OEIS entry, 
the initial terms of the sequence are correct. For n — 1, ... ,8 they are 



More generally, let kbe a positive integer not divisible by p, where q is a power of p. We 
consider the solutions of A k = I and let a n be the number of n x n solutions with coefficients 
in Fq. Now z k — 1 factors into a product of distinct irreducible polynomials 



Note that a n counts the homormorphisms from the cyclic group of order k to GL n (g). The 
generating function given here is a special case of the generating function given by Chigira, 
Takegahara, and Yoshida [3] for the sequence whose nth term is the number of homomor- 
phisms from a finite group G to GL n (g) under the assumption that the characteristic of F q 
does not divide the order of G. 

Now we consider some specific examples. Let k = 3 and let q be a power of 2. Then 




1, 4, 22, 316, 6976, 373024, 32252032, 6619979776. 



For q = 4 the sequence is A053856 and the initial terms for n = 1, . . . , 7 are 



1, 16, 316, 69616, 21999616, 74351051776, 374910580965376. 



z k -l =(j) 1 (z)<f) 2 {z)---(i) r {z). 



Let di = deg0j. Then the generating function for the a n is 




z 3 - 1 = z 3 + 1 = (z + \){z 2 + z + l) 
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is the irreducible factorization. Thus, d\ = 1 and d,2 = 2. The generating function is 

^ 7n \ >n ^™ / V >n 7m(<7 ) / 

n>l \m>0 / \m>0 / 

With q = 2 we get the sequence A053725 with initial terms for n = 1, . . . , 8 
1, 3, 57, 1233, 75393, 19109889, 6326835201, 6388287561729. 
With q = 4 we get the sequence A053857 with initial terms for n = 1, . . . , 5 

3, 63, 8739, 5790339, 25502129667. 
For the next example, let k — 8 and let q be a power of 3. Then 

J _ 1 = ( z _ 1)(2 + 1)0 3 (2)0 4 ( Z )0 5 (2), 

where the last three factors all have degree 2. Thus, r = 5, d\ = d,2 = 1, and — — — 2. 
The generating function is given by 

EMEfVfe^y. 

n>l 7n \m>0 7m / \m>0 lmW V 

For q = 3 this gives the sequence A053853 beginning for n — 1, ... ,5 

2, 32, 4448, 3816128, 26288771456. 
These are the sequences in the OEIS for small values of k and q. 



k 


q = 2 


q = 3 


q = 4 


2 


A053722 


A053846 


A053856 


3 


A053725 


A053847 


A053857 


4 


A053718 


A053848 


A053859 


5 


A053770 


A053849 


A053860 


6 


A053771 


A053851 


A053861 


7 


A053772 


A053852 


A053862 


8 


A053773 


A053853 


A053863 


9 


A053774 


A053854 




10 


A053775 


A053855 




11 


A053776 






12 


A053777 







1.12 Nilpotent matrices 

Fine and Herstein |3] proved that the number of nilpotent n x n matrices is g n ( n_1 ), and 
Gerstenhaber [7j simplified the proof. Recently Crabb has given an even more accessible 
proof [4^ . For q = 2 this is sequence A053763 with initial terms from n = 1 to n = 6 

1, 1, 4, 64, 4096, 1048576, 1073741824. 
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1.13 Cyclic matrices 

A matrix A is cyclic if there exists a vector v such that {v4 l i;|z = 0,1,2,...} spans the 
underlying vector space. (The term regular is also used. ) An equivalent description is 
that the minimal and characteristic polynomials of A are the same. Let a n be the number 
of cyclic matrices over F q . The generating function factors as 

n>l !n d>l 

and can be put into the form 

4jt» ?{?-*)) 

where v& is the number of irreducible, monic polynomials of degree d over F q . The generating 
function can be extracted from the proof of Theorem 1 in which we present in section 2. 
For q = 2 the sequence from a\ to a-j is 

2, 14, 412, 50832, 25517184, 51759986688, 422000664182784. 

Wall proved that the probability that an n x n matrix is cyclic has a limit 




Fulman [0] also has a proof. For q = 2 this limit is 0.7403 



1 + 



u 



1 1 



(u/qY 



"el 



1.14 Semi-simple matrices 

A matrix A is semi-simple if it diagonalizes over the algebraic closure of the base field. Let 
a n be the number of semi-simple n x n matrices over F q . Then the generating function has 
a factorization 

where r )j{q d ) = \GL j(q d )\. For q = 2 the sequence from a\ to aj is 

2, 10, 218, 25426, 11979362, 24071588290, 195647202043778. 




1.15 Separable matrices 

A matrix is separable if it is both cyclic and semi-simple, which is equivalent to having 
a characteristic polynomial that is square-free. Let a n be the number of separable n x n 
matrices over F q . Then the generating function factors 

n>l ,n d>l V H 
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This can also be factored as 



n>l 




u d {u d - 1) 
q d (q d - 1) 




For g = 2 the sequence from ai to 07 is 

2, 8, 160, 22272, 9744384, 20309999616, 165823024988160. 

The number of conjugacy classes of separable n x n matrices is q n — q n ~ x for n > 2 and 
q for n = 1. This is proved by Neumann and Praeger fH\ Lemma 3.2] by showing that the 
number of square-free monic polynomials of degree n is q n — q n ~ 1 . There is a natural bijection 
between the set of conjugacy classes of separable matrices and the set of square- free monic 
polynomials obtained by associating to each conjugacy class the characteristic polynomial 
of the matrices in that class. 

1.16 Conjugacy classes 

Let a n be the number of conjugacy classes of n x n matrices over F q . The generating function 
for this sequence is 



which is the reciprocal of the limiting probability that a matrix is invertible. See section 1.2. 
For q = 2 the sequence is A070933. The initial terms for n — 1, . . . , 10 are 



Let b n be the number of conjugacy classes in the general linear group GL n (q), the group 
ofnxn invertible matrices over F q . Then 




The number of conjugacy classes grows like q n . In fact, 




2, 6, 14, 34, 74, 166, 350, 746, 1546, 3206. 



For q = 3 the initial terms for n = 1, . . . , 10 are 



3, 12, 39, 129, 399, 1245, 3783, 11514, 34734, 104754. 




For q = 2 the sequence is A006951, which starts 



1, 3, 6, 14, 27, 60, 117, 246, 490, 1002. 
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For q = 3 the sequence is A006952, which starts 



2, 8, 24, 78, 232, 720, 2152, 6528, 19578, 58944. 

For q = 4, 5, 7 the sequences are A049314, A049315, and A049316. 
The number of conjugacy classes is asymptotic to q n : 

lim — = 1. 

n^oo q n 

Hence, in the limit the ratio of the number of conjugacy classes of invertible matrices to the 
number of conjugacy classes of all matrices is the same as the limiting probability that a 
matrix is invertible. That is, 



lim ^ = nfi--V 

n^coa n A± ^ q rj 



Sequence A070731 gives the size of the largest conjugacy class in GL n (2). Starting with 
n — 1 the initial terms are 

1, 3, 56, 3360, 833280, 959938560. 

The minimal order of the centralizers of elements in GL n (2) is given by the quotient of 7 n 
by the nth term in this sequence. The resulting sequence for n = 1, ... 10 is 

1,2,3,6,12,21,42,84,147,294. 

This sequence is A082877 in the OEIS. 

2 Selected Proofs 

2.1 The cycle index and generating functions 

In sections 1.6-1.13 we make heavy use of generating functions of the form 



n>l 



-U 

In 



where the sequence a n counts some class of n x n matrices. These generating functions come 
from the cycle index for matrices that was first defined by Kung |Hj and later extended 
by Stong ^2]- We follow Fulman's notation in [Hj. The cycle index for conjugation action 
of GL„(g) on M n (q) is a polynomial in the indeterminates x^ t \, where ranges over the 
set $ of monic irreducible polynomials over F q and A ranges over the partitions of the 
positive integers. First we recall that the conjugacy class of a matrix A is determined by the 
isomorphism type of the associated F q [z]-module on the vector space F™ in which the action 
of z is that of A. This module is isomorphic to a direct sum 

9©F t H/(^) 

i=\ j=i 
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where <pi, . . . , are distinct monic irreducible polynomials; for each i, A^i > Aj^ > • • • > Xi,^ 
is a partition of rij = ^ . Ajj. Since A is n x n, then n = ^ nj deg 0j = . Ajj deg 0,. Let 
Aj denote the partition of rii given by the Xij and define |Aj| = rii. The conjugacy class of A 
in M n (g) is determined by the data consisting of the finite list of distinct monic irreducible 
polynomials <f>i, . . . , (j>k and the partitions Ai, . . . , A&. 
The cycle index is defined to be 



where \<f>(A) is the partition associated to <ft i n the conjugacy class data for A. If <p does not 
occur in the polynomials associated to A, then \<f>(A) is the empty partition, and we define 

We construct the generating function for the cycle index 

«>! ,n AeM„(«)*£* 

This generating function has a formal factorization over the monic irreducible polynomials. 
For the proof we refer to the paper of Stong [T2J p. 170] . 

Lemma 1 

v) ri.n ™ „,|A|deg<A 

e n^,=nE^r- 

where c^(A) zs i/ie order of the group of module automorphisms of the F q [z]-module ® . F g [^]/ (0 
Lemma 2 Fzx a monic irreducible <p. Then 



A, 



7 r deg <f> 



a c ^ A ) 4i v g ' 

Proof We use the formula for c^(A) proved by Kung [HI Lemma 2, p. 146]. Let A be a 
partition of n and let hi be the number of parts of size i. Let di = b\ + 2&2 + • • • + z&j + 
i(fo i+ i H Vb n ). Then 

C^(A) = \\ JJ( g *dog0 _ g (d s -fe)deg^_ 
i fc=l 

We see from this formula that c^(A) is a function of A and g deg ^. Therefore it is sufficient to 
prove the lemma for a single polynomial of degree 1, say <ft(z) = z, and then to replace u by 
u &eg<p an j q gdeg<^ g Q ^ now we w jjj p rove that for (f)(z) = z, 



? c ^) Si 1 <T 
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We split the left side into an outer sum over n and an inner sum over the partitions of n 

y\\\ -.71 



To evaluate the inner sum we note that 7„ X}|A|=n ^Jxj * s the number of n x n nilpotent 
matrices, which is g ra(n_1 ) from the theorem of Fine and Herstein j^j. Therefore, the u n 
coefficient of the left side is 



y — 



l 



For the right side we use a formula of Euler, which is a special case of Cauchy's identity and 
a limiting case of the g-binomial theorem. A convenient reference is the book of Hardy and 
Wright (HI Theorem 349, p. 280]. For \a\ < 1, \y\ < 1, the coefficient of a n in the infinite 



product 



ay r ) 

r>l v y 1 



IS 



(l_ y )(l_ y 2)... (1 _ y n)- 

Let a = u and y = - to see that the u n coefficient of 



r>l 

is 



i 

u 



1 



Therefore, the coefficients of u n of the left and right sides are equal and so we have proved 
that for <p(z) = z, 



□ 

Lemma 3 Let $' 6e a subset of the irreducible monic polynomials. Let a n be the number of 
n x n matrices whose conjugacy class data involves only polynomials <fi e T/ien 

„ „,|A|deg</> 

n>l A rv 

a„ „ T-r T-r A M d6g0 



n>l m </>g$'r>l v ^ 
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Proof Using Lemma 1 we set x^x = 1 or according to whether G $' or not. Then on 
the left side the inner sum is over the matrices A that do not have factors of <fi G $' in their 
characteristic polynomials, and so the coefficient of u n is simply the number of such n x n 
matrices. This gives the first equality in the statement of the lemma 

. n . „, | A| deg 

n>l ln <j>^> A C ^ A) 

Then from Lemma 2 we get the second equality. □ 
Lemma 4 

u \\\dcg(f> 



nz 



nn 1 



l-u 11 11 V a rde ^ 

Proof Using Lemma 3 we let $' be the complement of the polynomial <f>(z) = z. The ma- 
trices not having factors of z in their characteristic polynomials are the invertible matraices. 
Thus, a n = 7„, and so the left side is 1 + J2 n >i M ™- n 
We can generalize the first part of Lemma 3 to allow for conjugacy class data in which the 
allowed partitions vary with the polynomials. The proof follows immediately from Lemma 
1. 

Lemma 5 For each monic, irreducible polynomial (f> let be a subset of all partitions of 
the positive integers. Let a n be the number of n x n matrices A such that X^A) G for all 
(p. Then 

. n . „ |A|dcg0 

n>l ln 4> ASL °^ A} 

Lemma 6 

,deg<p\ / d \ v d 



d> H 7 d>l V ^ 7 



Proof Unique factorization in the ring of polynomials F q [z] means that each of the q n 
monic polynomials of degree n has a unique factorization as a product of monic irreducible 
polynomials. This implies the factorization of the generating function 



E^riE 



yk deg < 



n>0 </> k>0 

The left side and the inner sum on the right are geometric series, and so 



— =n— 

1 — qu A A 1 — u 
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Grouping the factors on the right according to degree we get 

,d 



— =n(— 

l- q u 11 Vl- 
Substituing -u/g for u gives 

T^ = n(i-( M /g) d ) ' 

The lemma follows by taking reciprocals. □ 
The final lemma below is used to evaluate limiting probabilities such as \im n ^ 00 a n /'-y n . 
The proof is straightforward and omitted. 

Lemma 7 // 

z — ' 1 — u 

n>0 

where F(u) is analytic and the series for F '(1) is convergent, then 

lim a n = F(l). 

n^oo 

Now with these lemmas we are ready to prove several results stated in §1. 
2.2 Linear and Projective Derangements 

In this section we prove the results given in §1.8. Linear derangements are matrices with no 
eigenvalues of or 1, which means that their characteristic polynomials do not have factors 
of z or z — 1. 

Theorem 8 Let e n be the number of n x n linear derangements. Then 

n>l ln 1 U r>l \ q ' 

Proof In Lemma 3 let $' = $ \ {z, z — 1} to see that 

n>l ,n r>l \ ^ / 

On the right side multiply and divide by the product corresponding to <p(z) = z, which is 

r>l v ^ 

to see that 

HE^=nn('-3"n('-; 

n>l 7n «^ 2 r>l V y / r > x V y 
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Then use Lemma 4 to give 



n>l r>± 



U 
7 



From the generating function we can derive a recursive formula for e n . 
Corollary 9 

e n = e n ^(q n - l)?"" 1 + (-1) V ( "~ 1)/2 , e = 1. 
Proof From Theorem 8 it follows that e n /7„ is the sum of the u 1 coefficients of 

Eft 1 "«/?') 

r>l 

for i — 0, 1, . . . , n. Now the u % coefficient is 

(-i)* y — - — • 

\<r\<T2<---<Ti 

By induction one can easily show that this coefficient is 

f-iV 



Therefore 

7n ^ B (9*- IK?*" 1 -!)•••(? -1)' 

Next, 

e n _ e n _i (-l) n 



□ 



7n 7n-l (9 n -l)---(5-l)' 

Making use of the formula for 7„ and 7 n _i and canceling where possible we see that 

e n = e n ^(q n - l)^" 1 + (-1) V (n_1)/2 . 

□ 

We present the next result because the sequence e n for q = 2 is given in the OEIS as 
2 n ( n-1 )/ 2 a n where a„ satisfies a recursive formula. 

Corollary 10 Lei 

6 77. 



On. = 



*» q n(n-l)/2 ■ 

Then a n satisfies the recursion 

a n = a n . 1 (q n -l) + (-l) n , a = l. 
Proof The proof follows immediately from the recursive formula for e n . □ 
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Corollary 11 The asymptotic probability that an invertible matrix is a linear derangement 



IS 



The asymptotic that any matrix is a linear derangement is 

2 



"»-£=n(i-V) ■ 

n— >oo Q™ \ <7 / 

r>l x ' 



Proof We use Lemma 7 for the first statement. Then 

hm — j = l im 2 



(lim — ) ( lim ) . 
n-»oo ry n J yn^oo q n J 



We have just computed the first limit on the right and the second limit (from §1.2) is the 
same. □ 
The proof of Theorem 8 easily adapts to give the generating function for the number of 
projective derangements. 

Theorem 12 Let d n be the number of n x n projective derangements. Then 

9-1 



n>l ln 1 U r>l V 9 / 



Proof In this case $' = $ \ {z — a\a G F q }. As in the proof of Theorem 8 this time 
we multiply and divide by the products corresponding to all the linear polynomials except 
(j)(z) — z — 1. There are q — 1 of these and so we get 

n>l < n 4>+zr>\ \ y / r>l V ^ 7 

Use Lemma 4 to finish the proof. □ 
It would be interesting to find a recursive formula for the d n analagous to that given in 
Corollary 9 for the e n . 

2.3 Diagonalizable matrices 

Theorem 13 Let d n be the number of diagonalizable n x n matrices. Then 

d n „ fsr^ u rr X 



and 



n>0 7n \m>0 7m 



d n = 



In 



.+„„=„ 7ni " ' ln " 
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Proof We use Lemma 5. Diagonalizable matrices have conjugacy class data that only 
involves the linear polynomials <fi(z) = z — a, a e F q and partitions \<j,(A) that are either 
empty or have the form 1 > 1 > . . . > 1. These partitions are indexed by non- negative 
integer. For <f>(z) = z — a and for A the partition consisting of m l's, the corresponding 
mxm matrix in the canonical form is the diagonal matrix with a on the main diagonal. The 
centralizer subgroup of this matrix is the full general linear group GL m (g) and so c$(\) = j m . 
Then from Lemma 5 we see that 



It should be noted that the formula for d n follows directly from the knowledge of the 
centralizer subgroup of a diagonalizable matrix and does not require the full machinery of 
the cycle index generating function. Also, there is an alternative approach to this problem 
given in IT(T . 

2.4 Projections 

Theorem 14 Let p n be the number of n x n projections. Then 



Proof Since projections are diagonalizable matrices with eigenvalues restricted to be in the 
set {0, 1}, it follows that the product in the generating function for the cycle index is over 
the two polynomials z and z — 1 and the partitions are restricted as described in the previous 
theorem. Thus, we see that 




From this the formula for d n follows immediately. 



□ 




and 





The formula for p n follows immediately. 



□ 
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2.5 Solutions of A k = I 

Theorem 15 Let a n be the number ofnxn matrices A satisfying A 2 = I, and assume that 
the base field F q has characteristic two. Then 



a n = 



In 



g*(2n-3*) 7i 7 n _ 2i " 
0<i<n/2 

Proof The rational canonical form for A is a direct sum of companion matrices for z — 1 
and (z — l) 2 . No other polynomials occur. Thus, the partition X^(A) for <f>(z) = z — 1 consists 
of b\ repetitions of 1 and 62 repetitions of 2, where 61 + 26 2 = n. From Lemma 5 we see that 

dm - NT^ U bl+2b2 



where b = (61,62) denotes the partition. To compute c z _i(6) we refer to the formula stated 
in the proof of Lemma 2. In that formula we have d\ — b\ + 62 and d 2 = 61 + 26 2 , and then 

i=i fc=i 

This becomes ^ ^ 

c 2 -i(6) = f[(q dl ~ q dl - k ) f[(q d2 - 9*"*)- 
fc=i fc=i 

Now let i — d± — k in the first product and let i = d 2 — k in the second product, so that 

di-l d 2 -l 

c,-i(b) = II (q dl ~ ?*) II 

i=&2 i=d,2—b2 

From each factor of the first product we pull out a factor of q h2 and from each factor of the 
second product we remove a factor of q d2 ~ b2 . This gives 

c z -i{b) = q^ b2 lbl q^ d2 ' b ^ b2 lb2 . 

Then 

On _ v - ^ 1 

where we note that d 2 = n and 61 = n — 26 2 . Therefore, 

a n _ \ ^ 1 

^ ^h=n^ n ' 2h2) ^ b2 

which is equivalent to the formula to be proved. □ 
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Theorem 16 Assume that k is a positive integer relatively prime to q. Let 



z k -1 = 0i(*)02(*) •••&■(*) 



be the factorization of z k — 1 over F q into distinct irreducible polynomials with di = deg <pi 
and let a n be the number of n x n matrices over F q that are solutions of A k = I. Then 



r 



a 



•n 



n>l 



7n 



i=l m>0 



lm(q dl )' 



Proof The rational canonical form of a matrix A satisfying A k = I is a direct sum of any 
number of copies of the companion matrices of the 0j. Thus, with Lemma 5 the product is 
taken over the (pi for % — 1, . . . , r and the subset L^. is the same for all i and consists of the 
partitions in which all parts are 1. Hence, L^. can be identified with the non- negative integers 
m = 0, 1, 2, . . .. For the partition A given by m l's, the value of c^(A) is the order of the 
group of automorphisms of the Fjzj-module Fjz]/ (<fi™) , but this module is the direct sum of 
m copies of the extension field of degree di over F q . Therefore, the group of automorphisms 
is the general linear group GL m (q di ), and so c^X) = r y m (q di )- □ 

2.6 Cyclic matrices 

Theorem 17 Let a n be the number of cyclic n x n matrices. Then 



Proof In terms of the conjugacy class data a matrix is cyclic if has at most one part for 
each 0. Thus, L^ = {0, 1,2,.. .}, where m means the partition of m having just one part. 
For these partitions we have c<^(0) = 1 and c^(m) = g mdc s<^ — g(m-i)deg<?>_ ttn the 
number of cyclic matrices of size n x n. Beginning with Lemma 5 we find the generating 




and 
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function. 



111 o mdeg0 o (m-l)deg0 

<* \ m>l H H 



n i+E 



•md ' 



d>l \ m>l y y 



<2>1 \ y m>l H 



u d ^ 



For the second part, we use Lemma 6 to get 

1+ ' 



n>l ,n d>l v ^ 7 <i>l 



g d - 1 1 - (u/q) d / 

Combine the products and simplify to finish the proof. □ 
2.7 Semi-simple matrices 

Theorem 18 Let a n be the number of n x n semi-simple matrices over F q . Then 

( i,0d X Vi 

IS 7 " ii V 

Proof In order for a matrix to be semi-simple the partition associated to a polynomial 
has no parts greater than I. Thus, for all 0, L<p = {0, 1, l 2 , . . . , P, . . .}, where 1 J means the 
partition of j consisting of j copies of 1. Then c^V) is the order of the automorphism group 
of the F 9 [,z] -module which is the direct sum of j copies of F q [z]/((f>), which can be identified 
with (Fgd) J , where d = deg</>. Furthermore, the automorphisms of the module (F 9 [z]/(0)) ej 
can be identified with the F g d-linear automorphisms of (F q dy . This group of automorphisms 
is GL j(q d ), and so c^(P) = , ~fj{q d )- 



From Lemma 5 we find 

On , TT \ ^ II 

In 



I A I dcg r 



n>l 




n(i + E 7j(9 « 
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□ 



2.8 Separable matrices 

Theorem 19 Let a n be the number of separable n x n matrices over F q . Then 



n>l ln d>l V ^ 



and 

u d {u d -l)^ Ud 



n>l 7n 1 w d>l v 



- 1) 



Proof For a matrix to be separable the allowed partitions are either empty or the unique 
partition of 1. So the sum over in Lemma 5 is a sum of two terms. For A = 0, 



and for A = 1, 



Therefore, 



u \X\deg(p 

c^(A) 

u \X\dcg<f> u deg< ' 



C^(A) g*S*-l" 



„ „,|A|deg<; 



n>l 



gdeg</> _ ]_ 



„>i v « 

For the second statement of the theorem multiply the right side in the line above by 

d\ u d 



d>i v y 7 



which is equal to 1 by Lemma 6. Then combine the products. □ 
2.9 Conjugacy classes 

Theorem 20 Leta n be the number of conjugacy classes ofnxn matrices. Then the ordinary 
generating function for the sequence {a n } is given by 



1 + Ev"IIr^' 
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Proof Recall the ordinary power series generating function for partitions factors as the 
infinite product 

n>0 r>l 

where p n is the number of partitions of the integer n. A conjugacy class ofnxn matrices 
is uniquely specified by the choice of a partition for each monic irreducible polynomial 
such that J2s \ ^<t> \ deg0 = n. Consider the infinite product over <p 



radcg ( 

n>0 



The coefficient of u n is a sum of terms of the form p ni u ni deg ^ • • ■ p nk u nkdes ^ k where n = 
^nfcdeg^fc. Therefore, the ^"-coefficient is the number of conjugacy classes ofnxn ma- 
trices. Then 

n>l <f, n>0 



,r deg 4> 

r>l 4> 

\ Vd 

r>ld>l 

Then by substituting qu r for u in Lemma 6 and inverting we see that 



nn^ 

4, r>l 
r>l cf> 



d>i 

Therefore 



n a 



u rd J 1 — qu 
1 + V a n u n = TT ; 

n>l r>l * 



□ 
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